Jr 


n. 


— 


ROBERT  G.  BURETTE 
Major,  USAE 


Technical  Report  RK-CR-74-13 


2±Li 


ANALYSIS  OF  CHEMICAL  LASERS 

Volume  3 

One-Dimensional  Laser 
and  Mixing  Program  Guide 


v-  a 


sj 


m 

o 


D  ■=—  CD 

— —  5 

SS5  * 

3 


ru 


Larr/  R.  Ring 
Ronald  R.  Mikatarian 


LockheedMissU^^&rSpac^Cfcjijj^y^jyi^ 
T^untsvTn^ire^eS^Sn^lInginlering  Center 
4800  Bradford  Drive 
Huntsville,  Ala.  35807 


June  1974 


Final  Report  for  Contract  DAAI01-74-C-0173 


!  OAN  COPY:  RETURN  TO 
"RL/PSTL  TECHNICAL  LIBRARY 
TTLAND  AFB,  NM  87117-5776 


[  Appro 


Approved  for  public  release; 
\  distribution  unlimited 


Prepared  for 


20080303056 


Propulsion  Directorate,  U  S  Army  Missile  Research 
Development  and  Engineering  Laboratory 
U  S  Army  Command,  Redstone  Arsenal,  Alabama  35809 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as 
an  official  Department  of  the  Army  position  unless  so 
designated  by  other  authorized  documents. 


DISPOSITION  INSTRUCTIONS 

Destroy  this  report  when  it  is  no  longer  needed.  Do 
not  return  it  to  the  originator. 


TRADE  NAMES 

Use  of  trade  names  or  manufacturers  in  this  report  does 
not  constitute  an  official  endorsement  or  approval  of  the 
use  of  such  commercial  hardware  or  software. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 


TECH  LIBRARY  KAFB,  NM 


REPORT  DOCUMENTATION  PAGE 


1  REPORT  NUMBER 

RK-CR-74-  1  3 


2.  GOVT  ACCESSION  NO. 


READ 

BEFORE 


_  0147102 

3  RECIPIENT'S  C A  l  ALOG  NUMBER  I! 


Title  (nn<i  Subtitle) 

Analysis  of  Chemical  Lasers,  Vol.TII, 
One-Dimensional  Laser  and  Mixing  Program 
Guide 


5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Report 

Oct  1973  -  June  1974 


6.  PERFORMING  ORG.  REPORT  NUMBER 

LMSC-HREC  TR  D390222-III 


7.  AuTHORfsJ 


Larry  R .  Ring 
Ronald  R.  Mikatarian 


8.  CONTRACT  OR  GRANT  NUMBERfs) 


DAAHO  1-74-C-0173 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Lockheed  Missiles  &  Space  Company 
Huntsville  Research  &  Engineering  Center 
Huntsville,  Ala.  35807 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

Project  IT562609A308 
AMC  MS  662609-55-4730 


11-  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Propulsion  Directorate,  US  Army  Missile 
Research,  Development  and  Engineering 
Laboratory,  US  Army  Missile  Command, 
Redstone  Arsenal.  Ala.  35809 


12.  REPORT  DATE 

June  1974 


13.  NUMBER  OF  PAGES 


U  MONITORING  AGENCY  NAME  8  ADDR  ESSfi/  different  from  Controlling  Office) 

Same  as  (11)  above 


15.  SECURITY  CLASS,  (of  thla  report) 

Unclassified 


15«.  DECLASSIFICATION  'DOWNGRADING 

SCHEDULE 

N/A 


16  DISTRIBUTION  STATEMENT  (of  thla  Report) 

Approved  for  Public  Release;  Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 

Same  as  (16)  above 


18  SUPPLEMENTARY  notes 


N/A 


19  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  Identify  by  block  number) 

Chemical  Laser  Computer  Program 
Chemically  Reacting  Flows 
One-Dimensional  Flow 
Laser  Radiation 


20  ABSTRACT  f Continue  on  reverse  side  If  necessary  and  identify  by  block  number) 

The  One-Dimensional  Laser  and  Mixing  Program  (ODLAMP)  is  de¬ 
scribed  for  analyzing  multi- species  lasing  in  a  chemically  reacting  one¬ 
dimensional  stream  with  mass  addition  Mass  is  assumed  to  be  injected 
(at  a  predescribed  rate)  from  each  of  two  separate  streams,  into  a  primary 
lasing  stream.  Lasing  due  to  P-transitions  (i.  e.  ,  transitions  proceeding 
from  (v+  1 , J  -  1 )  to  (v,J))  is  treated.  Both  the  analysis  and  computer  pro¬ 
grams  are  presented,  including  an  input  guide,  sample  calculation  and 
flow  chart. 


DD 


FORM 
1  JAN  73 


1473 


EDITION  OF  1  NOV  85  IS  OBSOLETE 


SECURITY  CLASSIFICATION  of  THIS  PAGE  nihen  Data  Entered) 


FOREWORD 


This  document  is  Volume  III  of  a  five-part  final  report  which  presents 
the  results  of  work  performed  by  the  Lockheed  Huntsville  Research  &  Engi¬ 
neering  Center  under  Contract  DAAHO 1 -74 -C -0 173  for  the  Propulsion  Direc¬ 
torate,  U.  S.  Army  Missile  Research,  Development  and  Engineering  Labora¬ 
tory,  U.  S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama.  This  work 
was  monitored  by  Mr.  William  D.  Martin  of  the  Propulsion  Directorate. 

The  period  covered  by  this  report  was  from  20  September  1973  to 
30  June  1974. 

The  final  report  for  this  contract  is  comprised  of  the  following  volumes: 

1.  Laser  and  Mixing  Program  Theory  and  User's  Guide 
(LMSC-HREC  TR  D390222-I) 

2.  Chemical  Laser  Flow  Analysis  (LMSC-HREC 
TR  D390222-II) 

3.  One-Dimensional  Laser  and  Mixing  Program 
Guide  (LMSC-HREC  TR  D390222-III) 

4.  Method  of  Characteristics  Laser  and  Mixing  Program 
Theory  and  User's  Guide  (LMSC-HREC  TR  D390222-IV) 

5.  Rotational  Relaxation  Effects  (LMSC-HREC  TR  D390222-V). 

This  report  documents  the  work  done  to  modify  and  update  the  capabilities 
of  OD-LAMP  and  supersedes  report  Technical  Report  RK-CR-73-2.  A  summary 
of  the  modifications  incorporated  in  this  report  may  be  found  in  the  Introduction. 
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NOMENCLATURE 


Symbol 


A 

B 


C 


2 


f 

F. 

1 

AG 

h 

I 

K 

P 

kf 

£ 


reaction  rate  constant  (cm-particle- sec  units) 
activation  energy  (cal-mole  *) 

-2 

gravitational  constant  (980.6  cm-sec  ) 

-  8 

numerical  conversion  constant  (2.389  x  10 
cal-  sec -gm"  ^  -cm” 

catalytic  species  weighting  factor 

Y./W.  (mole.-gm"1) 

Gibbs  free  energy  (cal-mole  *) 
enthalpy  (cal-mole  *) 

/  -2  - 1  v 

internal  radiation  intensity  (cal-cm  -sec  ) 

equilibrium  constant 

reaction  rate  forward  rate  coefficient 
width  of  optical  cavity  (cm) 


L 


m 

P 

Q 

r 

o 

$ 

R 


T 

W 

u 


X 

Y. 

i 


width  of  stream  (cm) 

flow  rate  (gm-cm  *-sec  *) 
pressure  (torr ) 

radiation  contribution  to  energy  balance  (cal-cm”  *-sec”  *) 
mirror  reflectivities 

universal  gas  constant  (8312.97  dyne-cm-g  mole  *-K  *) 

species  production  (deletion)  due  to  stimulated  emission 
(mole.-cc"  1-  sec"  ^) 

time  (sec) 
temperature  (K) 
molecular  weight  (gm-mole”  * ) 
velocity  (cm-sec  *) 

species  chemical  production  (deletion)  rate 
(mole.-cc  -sec  *) 

longitudinal  distance  (cm) 
mass  fraction 
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NOMENCLATURE  (Continued) 


Greek 


a 

active  medium  gain  (cm 

_ 

threshold  gain  (cm  *) 

a 

P 

defined  by  Eq.  (10a) 

€ 

molar  photon  energy  (cal-mole  *) 

6m 

injection  flow  rate  (gm-cm  *-sec~ 

P 

density  (gm-cc’^) 

Subscript 

i 

r  ,  .th 

refers  to  i  species 

m 

refers  to  average  value 

L 

refers  to  lower  laser  level 

U 

refers  to  upper  laser  level 

1 

refers  to  initial  integration  station 

2 

refers  to  final  integration  station 

Superscript 

—  denotes  average 

I  refers  to  secondary  stream  (stream  2) 

II  refers  to  secondary  stream  (stream  3) 
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Section  1 
INTRODUCTION 

The  importance  of  continuous  wave  (cw)  chemically  pumped  lasers  has 
generated  a  number  of  investigations  designed  to  enhance  the  understanding 
of  laser  devices  and  their  potential.  These  investigations  have  led  to  the 
development  of  sophisticated  chemical  laser  models  (Refs.  1  through  3)  which 
couple  the  fluid  mechanic,  nonequilibrium  chemistry  and  laser  radiation  pro¬ 
cesses.  Unfortunately,  however,  these  computer  models  require  a  large 
amount  of  computer  time  to  analyze  one  set  of  conditions.  Consequently, 
when  performing  design  studies  where  many  variables  must  be  investigated 
the  computer  time  requirements  are  considerable. 

In  this  report  a  one-dimensional  chemical  laser  model  is  presented 
which  approximates  the  two-dimensional  effects  of  mixing.  Since  it  is  one¬ 
dimensional  and  requires  solution  along  only  one  streamline  as  opposed  to  the 
10  to  20  streamlines  required  in  two-dimensional  mixing  analyses,  computer 
run  times  are  reduced  significantly. 

As  Fig.  1  shows,  mass  from  two  dissimilar  secondary  streams  is  in¬ 
jected  into  a  primary  lasing  stream  at  some  prescribed  rate  which  varies 
with  longitudinal  distance.  Excitation  of  the  active  medium  traversing  a 
Fabry-Perot  optical  cavity  is  achieved  by  the  highly  exothermic  nonequilibrium 
reactions  occurring  in  the  cavity  region. 

The  program  can  handle  three  simultaneously  lasing  species  each  of 
which  in  turn  can  have  an  arbitrary  number  of  vibrational  levels.  Lasing  at 
either  fixed- J  or  shifting -J  rotational  lines  can  be  treated.  The  transition 
gains  (from  v+1,  J-l  — ►v,  J)  are  computed  accounting  for  the  effects  of  both 
Lorentz  and  Doppler  broadening.  In  addition,  the  program  can  handle  mirror 
absorptivities  and  reflectivities  which  vary  with  longitudinal  position. 
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Flow  Schematic 


The  program  allows  any  chemical  reaction  mechanism  to  be  investi¬ 
gated  as  long  as  the  rate  coefficients  and  species  thermodynamic  data  are 
known.  Thermodynamic  data  taken  directly  from  the  JANNAF  thermo¬ 
chemical  tables  is  input  in  tabular  form. 


In  this  report  the  governing  equations  and  a  complete  description  of 
the  computer  program  are  given  including  an  input  guide,  and  a  sample 
calculation.  The  output  of  the  computer  program  gives  the  longitudinal 
distribution  of  pressure,  velocity,  temperature,  species  mole  fractions, 
laser  transition  gains,  radiative  intensities,  and  laser  power  output. 


Modifications  made  and  incorporated  in  this  report  are  as  follows: 


•  Page  2-2  —  Equations  2a,  2b  and  2c  have  been  updated  to 
include  the  wall  for  ce  term  in  the  momentum  equation. 

This  term  was  inadvertently  not  included  in  the  original 
document. 

•  A  new  Section  7  —  Program  Usage  and  Comparison  has 
been  added.  This  section  briefly  discusses  the  program 
operational  experience  to  date  and  shows  a  comparison 
between  the  ODLAMP  and  LAMP  predictions  for  the  same 
sample  case. 

•  Page  B-5  —  The  input  has  been  modified  to  only  require 
input  of  the  collision  broadening  constants  (Card  20)  when 
Lorentz  and  Doppler  broadening  are  to  be  used  (LFLAG  =  1, 
Card  1). 

•  Appendix  C  —  Additional  sample  input  and  output  have  been 
provided  for  the  multi- specie  lasing  case  (DF  and  CO^). 
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Section  2 
ANALYSIS 


It  is  assumed  that  the  injection  rates  for  each  of  the  two  secondary 
streams  into  the  primary  lasing  stream  are  known  as  a  function  of  longi¬ 
tudinal  position.  The  computer  program  then  solves  for  the  fluid  mechanic 
and  chemical  properties.  The  flow  properties  (pressure,  velocity,  tempera¬ 
ture,  species  concentrations  and  cross  sectional  area)  are  defined  using  the 
conservation  of  momentum,  energy  and  mass  relationships  along  with  the 
continuity  equation  and  equation  of  state.  Closure  is  obtained  with  the  condi¬ 
tion  that  one  flow  parameter  (either  pressure,  temperature  or  cross-sectional 
area)  is  known. 

Referring  to  Fig.  1  the  governing  equations  for  the  three  streams  are: 


Mass 


Stream  (1)  —  Primary1' 


(p  u) 


rJ± 

A  ox 


w.  +  R. 

X2  X1 


Stream  (2)  —  Secondary 


F.  =  constant 
i 


Stream  (3)  —  Secondary 

I  I 

=  constant 


(la) 


(lb) 


(lc) 


In  Eq.  (la)  note  that  w  is  evaluated  at  the  forward  integration  point  (station  2). 
This  is  because  an  implicit  scheme  is  used  to  analyze  the  chemistry  effect. 
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Momentum 


Stream  (1) 

it 

^2  u2  =  ui  +  6m  ui  +  6m  ui  "  ci  (P2L2  'P1LP  +ClPm  ^L2"L1^  (2a) 


II  It 


where 
Stream  (2) 


D  =  P1  +  P2 


i  i  f 


™2  u2  =  ™1  u'l  "  6m'  U1  '  C1  (p2L2"  p'l  L'l)+  Clpm  (L2  "  Ll)  (2b) 


Stream  (3) 


m2  U2  =  ml  U1 


-  6m  Uj  -  C^p^L"  -  p''L»)  +ClPm(L2-  Lj)  (2c) 


Energy 


Stream  (1) 


1 


1 


m,  I  ■=—  +  C_  —=r  I  +  6m 


ll 


,W 


1 


2  2 


n 


1 1 2 


+  6m'  I  — +  C 


u. 


W, 


2  2 


1 


,  ,2 

h  u 

Jr-  +  c  2  -4- 

Wj 

=  m?  +  C, 

2\W2 


(3a) 


u 


~T  ]+  Q1(x2-x1)l 


Stream  (2) 


.2 

h  u 

4-  +  c2  ~Y~ 


W 


1 


4-+  C 

w. 


u. 


2  2 


(3b) 


Stream  (3) 


ii2' 


u. 


+  C 


,.2 


\W, 


2  2 


W. 


2  u2 
"  +  °2  2 


(3c) 
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Continuity- 


Stream  (1) 

.  .  .  ,  i  n 

m2  =  ml  +  + 

Stream  (2) 

m^  =  m^  -  6m 

Stream  (3) 

.  " 

m2 

where  m  =  puL 
State 

P 

W 


#  If  || 

=  -  6m 


Threshold  Condition 


a(x) 

Gain  (Along  Streamline) 


£n  (r  rT  ) 
o  L; 

2^ 


a.  T 
i=L 


€  Lp  (SU,L  FU  "  SL,  U  FL) 


(4a) 


(4b) 


(4c) 


(5) 


(6) 


(7) 


where  S^  L  and  S^  y  are  given  by  Eqs.  (2.29)  and  (2.30)  in  Ref.  1  for  the  case 
of  Doppler  and  Lorentz  broadening.  In  the  expressions  given  above  the  non- 
superscripted  parameters  refer  to  properties  in  the  primary  stream.  The 
constants  Cj  and  have  been  included  in  order  to  make  the  equations  dimen¬ 
sionally  consistent  based  upon  the  units  given  in  the  Nomenclature.  The  param¬ 
eters  R  and  Q  account  for  lasing  influences  upon  the  mass  and  energy  relations, 
respectively,  and  are  equal  to  zero  when  a.  <  a  .  The  subscripts  U  and  L  refer 
to  the  upper  and  lower  laser  levels,  respectively. 
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The  species  mass  fractions  computed  from  Eq.  (la)  do  not  take  into 

account  the  mass  injected  between  station  Xj  and  x^,  consequently  they  must 

be  modified.  This  is  done  as  follows:  The  mass  fractions,  Y.,  are  first  com- 

1 

puted  at  station  x^  from  Eq.  (la).  Then  they  are  modified  to  account  for  the 
mass  injected  between  Xj  and  x^  from 


m. 

1 


m  j  +  (6m  Ax) 

j  =  l 


where  Y'.  is  the  mass  fraction  computed  using  Eq.  (la). 
*2 


(8) 


The  parameters  R  and  Q  are  functions  of  the  gain,  a,  and  internal 
intensity,  I,  and  are  given  by 


and 


where 


R. 


“L1^  aL-lIL-l 


i=L  eL  eL- 1 


LS  VL 

°=  y  y  ,i. , 

L-t  J,  k  j,  k  j ,  k 

j  =  l  k=l 


h.  .  +  c .  ,  -  h  .  ,  ,  , 

p  =  j,k  j,k  j,k+l 


j»k 


Cj»  k 


(9) 


(10) 


(10a) 


th 

The  double  subscript  denotes  the  j  lasing  species  and  its  k  vibrational 
level  where  LS  is  the  number  of  lasing  species  and  VL  the  number  of  vib¬ 
rational  levels  for  each  lasing  species.  The  internal  intensities  are  obtained 
by  solving  a  set  of  simultaneous  equations  of  the  form 
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lh-l  GL-1  '  *L  GL  +  ^+1  GL+1 


VL 

E  *kG 

k  =  l 


daL 

LK  =  “cbT  '  Go 


(11) 


where  the  Gs  are  given  by  Eqs.  (3  5)  through  (39)  in  Ref.  1.  A  separate  set  of 
simultaneous  equations  must  be  solved  for  each  lasing  species. 

The  cavity  power  (for  each  lasing  line)  is  then  obtained  from 


P  .  (x)  =  /  I  z  dx 

out.  T  '  '  J  out. 
i=L  i 

xth. 

l 


(12a) 


where  (case  1) 


z  =  1 


(12b) 


when  the  optical  path  lies  in  the  plane  containing  the  three  streams, 
and  (case  2) 

L  +  L’  +  L" 


z  = 


(L  +  L-  +L")x=q 


(12c) 


when  the  optical  path  is  perpendicular  to  the  plane  containing  the  three  streams. 
In  Eq.  (12a) 


I  .  =  a.  I.Jt 

out.  T  11 
i=L 


(13) 


The  total  output  power(per  initial  width  of  three  streams  for  case  1  and  per  unit 
height  of  nozzle  bank  for  case  2)  is  then 


out 


LS  VL 

£  £  p 

j=l  i=l 


out. 


(14) 
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Section  3 

CHEMICAL  REACTION  RATE  EQUATIONS 
3.1  REACTION  RATES 


The  production  rates  for  all  chemical  species  are  calculated  in  the  usual 
fashion.  The  same  treatment  is  applied  to  the  various  vibrational  levels  of 
leacting  molecules,  and  of  the  lasing  molecules  in  particular,  that  is,  vibra¬ 
tional  levels  are  treated  as  individual  chemical  species. 

Twelve  types  of  reaction  or  vibrational  energy  transfer  mechanisms  are 
considered  as  possible  contributors  to  the  calculation  of  the  net  rate  of  produc¬ 
tion,  w.: 

1 


Reaction  type 


(1,  7) 

A  +  B 

(2,8) 

A  +  B  +  M 

(3,9) 

A  +  B 

(4,  10) 

A  +  B 

(5,  11) 

A  +  M 

(6,  12 

A  +  M 

— » 

C  +  D 

[1] 

- - t 

C  +  M 

M 

C  +  D  +  E 

[3] 

- - * 

C 

[4] 

~ - * 

C  +  D  +  M 

[5] 

- - * 

C  +  M 

[6] 

Reaction  types  (7)  through  (12)  correspond  to  reaction  types  (1)  through 
(6),  but  proceed  in  the  forward  direction  only. 


The  net  rate  of  production  for  all  reactions  is  given  below  in  the  form 
RP  ^  -  RM^  which  are  the  symbols  used  in  the  computer  program. 


1. 


kf  fa  fb 


k,  n2  F 
f  M  C  D 

K 

P 


(15a) 
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=  kf  p3  fafb 

kfp  fcfm 

M  K  ST 

P 

(15b) 

*«> 

“  V2fafb 

VJfcfdfe*t 

K 

p 

(15c) 

w«> 

■  kfp  fafb 

kf  P  Fc 

Kp  SRT 

(15d) 

w<J> 

=  k£p2FAFM 

kfp3FCFDFM5iT 

KP 

(15e) 

w«> 

=  kip2FAFM 

kfp2FCFM 

K 

D 

(15f) 

To  reduce  round-off  and  truncation  errors,  RP^  and  RM^  for  each 
reaction  are  computed  separately.  All  contributions  to  the  molar  rate  of  pro¬ 
duction  of  a  given  species  are  then  computed  and  added  algebraically  to  form 
wi*  Since  reaction  types  (7)through  (12)proceed  in  the  forward  direction  only, 
the  second  term  on  the  right-hand  sides  of  Eqs.  (15a)  through  (15f)  is  dis¬ 
regarded  in  calculating  the  contributions  to  w^. 

In  reactions  (2),  (5)and  (6)as  well  as  in  (8),  (11)  and  (12),  M  denotes  the 
catalytic  species  which  can  be  specified.  In  case  of  reactions  (2,  8),  (5,  11)  and 
(6,  12)the  situation  often  occurs  where  for  various  catalytic  species  the  re¬ 
spective  rate  constants  differ  only  by  a  constant  multiplier.  These  multipliers 
can  be  considered  as  third  body  efficiencies  or  weighting  factors.  If  such  a 

case  is  encountered,  the  third  body  species  mole  mass  ratio  F„  ,  becomes 

M 

effectively  a  fictitious  mole  mass  ratio,  consisting  of  the  weighted  sum  over 
all  those  species  having  a  nonzero  weighting  factor,  i.e., 


FM  "  £  fi  FM. 

i  i 


where  T  are  the  weighting  factors. 


(16) 


3-2 


3.2  RATE  CONSTANTS 


The  forward  rate  constant  is  generally  expressed  in  Arrhenius  form. 
The  equilibrium  constant,  K^,  is  determined  from  the  Gibbs  free  energy 
difference 


In  K  =  -  AG/KT 
P 


(17) 


For  speed  in  computation  the  rate  constants  are  divided  into  five  types: 


Rate  Coefficient  Type 


(1) 

< 

n 

'-W 

(18a) 

(2) 

-N 

kf  =  AT 

(18b) 

(3) 

kf  =  A  exp  (B/D?  T) 

(18c) 

(4) 

kf  =  AT"Nexp(B/!RT) 

(18d) 

(5) 

kf  =  AT  ~N  exp(B/»  TM) 

(18e) 
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Section  4 

THERMODYNAMIC  PROPERTIES 

To  make  maximum  usage  of  the  current  technology  of  finite  rate  formu¬ 
lations,  the  vibrational  levels  of  respective  excited  molecules  are  treated  as 
individual  chemical  species.  This,  of  course,  requires  separate  specification 
of  the  thermodynamic  properties  for  molecules  in  each  vibrational  level.  This 
is  accomplished  by  the  assumption  that  the  molecular  rotation  is  in  equilibrium 
at  the  translational  temperature  but  allows  the  molecules  to  vibrate  independently 
of  the  translational  temperature.  Accordingly,  each  molecule  is  permitted  to 
vibrate  independently  corresponding  to  the  energy  in  its  respective  vibrational 
level . 

Thermodynamic  data  is  input  directly  in  tabular  form.  The  thermo¬ 
dynamic  properties  (specific  heat,  entropy  and  fenthalpy)  for  each  species 
are  taken  directly  from  the  JANNAF  thermochemical  tables  (Ref.  4)  (or  other 
such  source)  and  input  to  the  program  as  described  in  Appendix  B.  Properties 
are  obtained  from  the  tables  by  linear  interpolations  on  temperature.  The 
Gibbs'  free  energy  is  computed  from  the  JANNAF  enthalpy,  entropy  and 
temperature . 
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Section  5 

PROGRAM  OPTIONS 


The  computer  code  has  been  developed  around  three  primary  options. 
They  are 


1.  The  longitudinal  distribution  of  temperature 

in  each  of  the  three  streams  is  known  a  priori. 

2.  The  longitudinal  distribution  of  pressure  in  each 
of  the  three  streams  is  known  a  priori, 

3.  The  longitudinal  distribution  of  the  total  cross- 
sectional  area  of  the  three  streams  is  known 

a  priori. 


Temperature  Option 


For  the  case  where  the  temperature  distribution  is  assumed  to  be  known 
the  calculational  procedure  is  as  follows: 


1.  Species  mass  fractions  are  computed  using  Eqs.  (1) 
and  (8). 

2.  The  velocities  in  the  primary  and  secondary  streams 
are  then  obtained  using  Eqs.  (3). 

3.  With  the  velocities  now  known  the  pressures  in  the 
three  streams  can  be  computed  using  Eqs.  (2). 

4.  The  densities  are  then  found  from  the  equations  of 
state  (Eq.  (5)). 

5.  The  areas  of  the  three  streams  are  computed  from 
continuity  (Eqs.  (4). 


Pressure  Option 


For  the  case  when  the  pressure  distribution  is  assumed  to  be  known 
the  calculational  procedure  is  as  shown  on  the  following  page. 
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1.  Same  as  Temperature  Option 

2.  The  velocities  are  computed  using  Eqs.  (2) 

3.  Temperatures  are  then  computed  using  the  energy 
equations  (Eqs.  (3)). 

4.  Same  as  Temperature  Option 

5.  Same  as  Temperature  Option. 


Area  Option 

When  the  area  is  assumed  to  be  known  the  solution  becomes 
iterative.  The  solution  follows  that  for  the  pressure  option  case.  A  pressure 
at  point  is  assumed  and  the  area  computed  and  checked  against  the  known 
area.  If  the  computed  and  known  areas  do  not  agree  the  assumed  pressure  at 
point  is  modified.  This  procedure  continues  until  agreement  is  reached. 

For  this  option  the  pressure  variations  for  all  three  streams  are  assumed 
to  be  identical,  but  unknown.  If  the  pressure  variations  are  not  taken  to  be 
identical  a  unique  solution  cannot  be  found. 
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Section  6 

PROGRAM  FEATURES 


•  Dynamic  Dimensioning 

In  order  to  establish  core  storage  requirements  which  are  sufficient  to 
run  a  particular  case  but  yet  keep  from  setting  aside  core  storage  which  is 
unused,  dynamic  dimensioning  has  been  incorporated  into  the  ODLAMP  code. 
Prior  to  compilation  of  the  ODLAMP  code,  therefore,  the  variable  subscripts 
used  in  the  common  and  dimension  statements  must  first  be  replaced  by  their 
appropriate  integer  values  for  the  particular  case  being  considered.  This  is 
done  by  running  the  ODLAMP  Dynamic  Dimensioning  Program.  This  program 
reads  as  input  the  ODLAMP  code  and  the  variable  subscript  names  and  their 
integer  values.  The  ODLAMP  code  is  then  searched  and  all  variable  sub¬ 
scripts  are  replaced  by  their  respective  integer  values,  thereby  generating  a 
new  ODLAMP  program  tape.  This  tape  is  then  rewound,  compiled  and  exe¬ 
cuted.  The  input  guide  for  the  ODLAMP  Dynamic  Dimensioning  Program  is 
given  in  Appendix  A. 

•  Plotting  (Page  or  4020  plots) 

If  so  desired  the  user  may  have  plotted  selected  parameters  (e.g.,  laser 
power,  pressure,  species  mole  fractions,  etc.)  as  a  function  of  distance.  This 
is  done  by  setting  IPLOT=0,  1,2  or  3.  If  IPLOT=0,  plots  are  omitted.  If 
IPLOT  =  1  page  plots  are  obtained  which  gives,  after  each  50  output  stations  as 
well  as  at  laser  cutoff,  a  plot  of  each  of  the  parameters  as  a  function  of  distance. 
^  £PLOT  =3  the  desired  data  will  be  stored  on  logical  unit  4  for  subsequent 
plotting  on  the  4020.  If  IPLOT  =2,  both  page  and  4020  plots  are  obtained.  The 

number  and  identification  of  the  parameters  to  be  plotted  are  read  in  on  cards 
22  and  23. 

•  Optical  Cavity  Width 

The  user  may  select  as  the  width  of  the  optical  cavity  either  the  width 
of  the  three  streams  which  is  evaluated  locally  or  a  width  which  is  input. 
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The  selection  of  which  option  desired  is  made  on  card  2.  If  the  width  is  input 
it  is  given  by  an  equation  of  the  form 

=  Q  +  Q.x  +  Q7x^ 
o  1  2 

where  and  x  are  in  centimeters.  The  coefficients  are  read  in  on  card  24. 
•  Mirror  Properties 

In  order  to  make  the  program  flexible, mirror  absorptivities  and  reflec¬ 
tivities  are  treated  as  variable  functions  of  distance.  They  are  given  by  equa¬ 
tions  of  the  form 


r  (or  a)  =  Qq  +  Q  +  (3  2x2  +  Q3x3  +  C^x"1  +Q5x"2  +  Q&  exp  (-x/Q?) 

The  coefficients  are  read  in  on  cards  8-11. 

•  Doppler/Lorentz  Broadening 

An  option  which  allows  the  user  to  select  either  Doppler/Lorentz 
broadening  or  Doppler  (only)  broadening  is  provided  (See  card  2).  However, 
it  is  recommended  that  the  Doppler/Lorentz  option  be  exercised  until  such 
time  as  the  user  becomes  familiar  with  the  conditions  where  the  Doppler 
(only)  broadening  approximation  is  adequate. 

•  Chemical  Kinetic  Rates 

An  output  option  is  provided  (card  2)  which,  if  exercised,  prints 
out  at  all  output  stations  the  species  kinetic  rates,  the  forward  and  reverse 
reaction  rates  for  all  chemical  reactions  and  the  forward  and  reverse  rate 
coefficients  for  all  chemical  reactions. 


6-2 


•  Mixing  Rate 


The  rate  at  which  mass  from  each  of  the  two  secondary  streams  is 
injected  into  the  primary  lasing  stream  is  assumed  to  be  known  a  priori  as 
a  function  of  longitudinal  position.  The  mass  injection  rates  are  then  de¬ 
scribed  by  an  equation  of  the  form 

6m  =  +  Q +  Q ^x2  +  Q^x3  +  Q4x"  1  +  Q5x“2  +  exp 

where  6m  and  x  are  in  gm-cm  *-sec  *  and  cm.  The  coefficients  are  read 
in  on  cards  6  and  7. 
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Section  7 

PROGRAM  USAGE  AND  COMPARISON 
7.1  PROGRAM  USAGE 

This  section  briefly  describes  some  of  the  operational  characteristics 
of  the  ODLAMP  program  which  have  been  observed  during  development  and 
check-out.  Relatively  little  operational  experience  has  been  obtained  to  date 
other  than  exercising  the  various  options  for  selected  test  cases.  Therefore, 
the  comments  that  follow  will  probably  not  cover  specific  characteristics 
which  will  be  uncovered  as  the  program  is  more  extensively  used. 

The  various  types  of  cases  considered  thus  far  may  be  summarized 
as  follows: 

•  Single  stream,  finite  rate  chemistry  analyses  in  variable 
area  ducts  without  laser  radiation. 

•  Simulation  of  multi -stream  viscous  mixing  of  dissimilar 
streams  with  finite  rate  chemistry  and  no  laser  radiation. 

•  Premixed  single  stream  analyses  of  one  lasing  specie 

(HF  or  CO,,)  with  both  fixed  and  shifting  J  levels  (Appendix  C-l). 

•  Premixed  single  stream  analyses  of  two  simultaneously 

lasing  species  (DF  and  CO2)  with  shifting  J  levels  (Appendix  C-2). 

•  Multi-stream  simulation  of  viscous  mixing  of  dissimilar 
streams  with  finite  rate  chemistry  and  one  lasing  species 
(HF). 

Program  execution  times  are  naturally  a  function  of  linear  step  size. 
Logically,  larger  input  step  sizes  would  mean  less  execution  time.  However, 
one  must  remember  that  the  laser  radiation  package  contained  in  the  program 
traces  discontinuities  such  as  threshold,  J  shift  and  cutoff.  The  program 
iteratively  solves  for  the  axial  location  of  these  discontinuities  with  a  high 
degree  of  accuracy.  The  actual  step  size  may  be  drastically  less  than  that 
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which  is  input.  Therefore,  specific  runtime  characteristics  are  difficult 
to  ascertain. 

Generally,  the  run  times  for  a  radiation  calculation  are  equivalent  to 
those  of  the  LAMP  code  for  a  premixed  stream  where  explicit  chemistry  is 
used.  However,  run  times  for  non-radiating  cases  are  drastically  reduced 
due  to  the  utilization  of  the  implicit  finite  rate  chemistry  model. 

Another  observed  characteristic  is  the  ability  of  the  program  to  handle 
the  transition  from  supersonic  to  subsonic  flow  conditions.  This  provides 
additional  utility  where  duct  flows  are  of  interest. 

Generally,  the  program  is  relatively  free  of  external  "fine  tuning  adjust¬ 
ments"  associated  with  many  numerical  techniques  and  requires  few  pre- 
execution  calculations  to  establish  input  data. 

7.2  ODLAMP  COMPARISON  WITH  LAMP  (REF.  1) 

The  ODLAMP  (one -dimensional  laser  and  mixing  program)  is  the  mani¬ 
festation  of  a  one -dimensional  theoretical  simplification  of  the  two-dimensional 
LAMP  (Ref.  1)  code.  The  ODLAMP  simulates  the  two-dimensional  viscous 
mixing  by  injecting  mass  at  prescribed  rates  into  a  primary  (lasing)  medium. 
Excitation  of  the  active  medium  traversing  a  Fabry-Perot  optical  cavity  is 
achieved  by  the  highly  exothermic  nonequilibrium  reactions  occurring  in  the 
cavity  region. 

A  sample  case  was  selected  and  the  data  input  to  both  the  LAMP  and 
ODLAMP  codes.  The  sample  case  was  selected  such  that  direct  comparisons 
of  the  flow,  chemistry  and  radiation  parameters  could  be  obtained.  A  hydrogen 
fluoride  (HF)  with  helium  (He)  diluent  medium  was  selected.  Six  vibrational 
levels  of  HF  were  considered  in  the  model.  The  following  table  contains  the 
input  initial  conditions. 
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U.  =  5.58  x  10^  m/sec 

P.  =  20  torr 

1  o 

T.  =  500  K 


Chemical  Species  Mole  Fraction 


HF(0) 


— HF(6) 

0.0 

F 

0.1 

F2 

0.01 

H 

0.0 

H2 

0.39 

He 

0.5 

Both  analyses  assumed  constant  pressure  as  a  function  of  axial  distance 
and  both  operated  at  the  rotational  level  producing  maximum  gain  via  J  shift¬ 
ing.  Fifty  reactions  and  five  catalytic  species  were  used. 

Figures  7-1,  7-2,  and  7-3  show  the  results  of  both  analyses.  Compari¬ 
sons  are  shown  of  the  axial  distributions  of  temperature  (Fig.  7-1),  species 
mole  fractions  (Fig.  7-2),  and  radiation  intensity  (Fig.  7-3)  for  both  the  ODLAMP 
and  LAMP.  This  check  case  demonstrates  the  ability  of  ODLAMP  to  duplicate 
the  LAMP  one -dimensional  premixed  calculations. 

Figure  7-4  illustrates  one  of  the  unique  features  of  ODLAMP,  i.e.  that 
of  simultaneous  lasing  from  two  different  molecules.  In  this  case  4  transitions 
of  DF  and  the  10.6 fj.  line  of  the  CC>2  molecule  were  all  handled  simultaneously. 

The  advantages  of  dynamic  dimensioning  used  in  ODLAMP  were  demon¬ 
strated  by  the  fact  that  ODLAMP  required  computer  core  storage  of  only 
57,0  00  octal  words  compared  to  134,000  octal  words  for  the  LAMP  code.  This 
becomes  a  definite  advantage  where  computer  systems  use  core  storage  re¬ 
quirements  in  assigning  job  priorities. 
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Fig.  7-1  -  Comparison  of  ODLAMP  and  LAMP  Predictions  of  Axial  Temperature  Distribution 


Specie  Mole  Fractions 


Fig.  7-2-  Comparison  of  ODLAMP  and  LAMP  Predictions  of  Axial  Species 
Mole  Fraction  Distribution 
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Internal  Intensity,  W/c: 


Fig.  7-3  -  Comparison  of  ODLAMP  and  LAMP  Predictions  of  Axial 
Radiation  Intensity  Distributions 
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Power  (kW/cm) 


Fig.  7-4  -  Calculated  Results  of  Simultaneous  Lasing  from 
Four  Transitions  of  DF  and  the  CO^  Upper  to 
Lower  Laser  Level 
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Appendix  A 

ODLAMP  DYNAMIC  DIMENSIONING  PROGRAM 
(Input  Guide) 


Appendix  A 


This  appendix  contains  an  input  guide  for  the  ODLAMP  Dynamic  Dimen¬ 
sioning  Program.  This  program  establishes  the  necessary  storage  require¬ 
ments  for  the  ODLAMP  code. 

There  are  seven  (7)  required  input  variables.  No  specific  order  of  input 
is  required.  The  data  are  input  in  the  format  of  one  variable  name  and  its 
corresponding  integer  value  per  data  card. 


Card 

Format 

Column 

Variable 

1 

110 

10 

NOP 

Number  of  input  variables 
(7  for  ODLAMP) 

2-8 

A10 

1-10 

OPT(l) 

Variable  name  (left  adjusted) 

110 

11-20 

OPT(2) 

Integer  value  assigned  to  the 
variable  name  (right  adjusted) 

The  following  is  a  list  of  the  required  input  variables  for  the  ODLAMP 
program. 

1  -  NS  =  Number  of  chemical  species  +1 

2  -  NM  =  Number  of  catalytic  species 

3  -  NMS  =  NS  +  NM 

4  -  NT  =  Number  of  temperature  values  in  thermo¬ 

dynamic  tables 

5  -  NR  =  Number  of  chemical  reactions 

6  -  NX  =  Total  number  of  lasing  transitions  for  all 

lasing  species 

=  NVS(l)  for  1  lasing  species 

=  NVS(l)  +  NVS(2)  +  1  for  2  lasing  species 

=  NVS(l)  +  NVS(2)  +  NVS(3)  +  2  for  3  lasing  species 

7  -  NY  =  NX  +  1. 
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Appendix  B 

ONE- DIMENSIONAL  LASER  AND  MIXING  PROGRAM 

(ODLAMP) 

(Input  Guide) 


Appendix  B 


This  appendix  contains  an  input  guide  for  the  One- Dimensional  Laser  and 
Mixing  Program  (ODLAMP)  developed  by  the  Lockheed-Huntsville  Research 
Engineering  Center.  This  guide,  along  with  the  information  from  the  rest  of 
this  report,  provides  the  necessary  information  for  the  user  of  this  program. 


Card 

1 

2 


Format 

Column 

Variable 

12A6 

1-72 

HDR 

1615 

5 

ITYPE 

9-10 

NT 

14-15 

19-20 

25 

NS 

NR 

KPWR 

30 

34-35 

39-40 

44-45 

50 

NM 

NV(1) 

NV  (2) 

NV  (3) 
LFLAG 

55 

IOUT 

60 

NSTRM 

65 

IGAS 

70 

KFLAG 

75 

IPLOT 

80  ILASL 


Description 


Title  card 

=  1  -  T  =  f(x)  known 
=  2  -  P  =  f(x)  known 
=  3  -  A  =  f(x)  known 
Number  of  temperature  points  in 
thermodynamic  tables 
Number  of  chemical  species 
Number  of  chemical  reactions 
=  0  Power  off 
=  1  Power  on 

Number  of  catalytic  species 

Number  of  lasing  transitions  for 
each  lasing  species 

=  0  -  Doppler  broadening 
=  1  -  Lorentz  +  Doppler  broadening 
=  0  -  No  species  +  reaction  rate 
printout 

=  1  -  Species  +  reaction  rate  print¬ 
out  (including  rate  coefficients) 
Number  of  streams  (3  max.  ) 

=  0  -  Gas  dynamics  and  chemistry 
only 

=  1  -  Includes  radiation 
=  1  -  Condensed  laser  output 
=  2  -  Detailed  output  for  each 
transition 

=  0  -  Omit  plots.  Do  not  read 
cards  22  and  23. 

=  1  -  Page  plots  only 
=  2  -  Page  and  4020  plots 
=  3  -  4020  plots  only 
=  0  -  Width  of  lasing  zone  based 
on  overall  width  of  streams 
=  1  -  Width  of  lasing  zone  read  in 
as  input.  Read  card  24. 


B-  1 


NOTE:  Cards  3,  4  and  5  contain  the  independent  variable  coefficients.  The 

independent  variable  type  is  defined  by  ITYPE  (Card  1).  The  physical 
units  input  are: 


ITYPE 

Independent  Variable 

Units 

1 

Temperature 

K 

2 

Pressure 

tor  r 

3 

Area 

cm 

The  known  axial  distributions  of  the  variables  on  Cards  3  through  11 
are  input  through  equations  of  the  form: 


Q  =Q0  +  Q 

jX  +  Q2x2  +  Q 

3x2  +  Q4x 

1  +  Q5x  2  +  exp  (_x/Q7) 

Card  Format 

Column 

Variable 

Description 

3 

8E10.6 

1-10,  etc. 

COEF(l) 

STRM1  (T,P,  A)  coefficients  (center)  ' 

4 

8E10.6 

1-10,  etc. 

COEF(2) 

STRM2  (T,P,  A)  coefficients  (upper) 

(Read  only  if  NSTRM  >  1) 

5 

8E10.6 

1-10,  etc. 

COEF(3) 

STRM3  (T,P,  A)  coefficients  (lower) 

(Read  only  if  NSTRM  =  3) 

6 

8E10.6 

1-10,  etc. 

COEF(4) 

STRM2  Mass  ejection  rate  (g/sec-cm) 
(Read  only  if  NSTRM  >1) 

7 

8E10.6 

1-10,  etc. 

COEF(5) 

STRM3  Mass  ejection  rate  (g/sec-cm) 
(Read  only  if  NSTRM  =  3) 

Input  Cards  8, 

9,  10  and 

1 1  only  if  IGAS  -1 - 

8 

8E10.6 

1-10,  etc. 

COEF(6) 

Mirror  1  reflectivity 

9 

8E10.6 

1-10,  etc. 

COEF(7) 

Mirror  2  reflectivity 

10 

8E10.6 

1-10,  etc. 

COEF(8) 

Mirror  1  absorptivity 

11 

8E10.6 

1-10,  etc. 

COEF(9) 

Mirror  2  absorptivity 

The 

following  cards  contain  the  thermodynamic  data*.*  The  first  card  contains 

the  species  name,  molecular  weight  and  heat  of  formation.  The  second  and 
remaining  cards  contain  the  temperature  and  corresponding  specific  heat, 
entropy  and  enthalpy  for  that  species.  Two  temperatures  and  corresponding 
thermodynamic  data  are  placed  on  each  card.  The  input  table  can  contain  up 
to  a  maximum  of  30  temperature  points.  The  data  are  input  exactly  as  pre¬ 
sented  in  the  JANNAF  tables. 


Stream  1  is  the  reacting/lasing  stream. 

The  order  of  the  species  must  be  identical  to  the  order  on  card  type  21a. 
Lasing  species  data  must  be  given  first,  in  the  order  v=0,  v  =  l,  v=2,  etc., 
followed  by  v=0,  v=l,  etc.,  for  the  second  lasing  species  (if  necessary)  and 
similarly  for  the  third  lasing  species. 
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Card 

F  ormat 

Column 

Variable 

Description 

12 

A6 

1-6 

AID 

Name  of  first  species 

El  0.3 

7-16 

WTMOLE 

Molecular  weight 

E10.3 

17-26 

HF 

Heat  of  formation,  1*293 .  (kcal/ 
mole)  1 

13 

E10.3 

1-10 

TTB 

* 

First  temperature  point  (K) 

E10.3 

E10.3 

11-20 

21-30 

CPTB 

GTB 

c_  (cal/mole-K) 

Sj  *  (cal/mole-K) 

E10.3 

31-40 

HTB 

hi  -  h298j  (kcal/mole)  ^ 

E10.3 

41-50 

TTB 

Second  temperature  point  (K) 

E10.3 

51-60 

CPTB 

Cp.  (cal/mole-K) 

Sj  (cal/mole-K) 

E10.3 

61-70 

GTB 

E10.3 

71-80 

HTB 

hj  -  1*293  (kcal/mole) 

NOTE:  1.  There  are  NS  card  no.  12s  and 

2*  NT/2  card  no.  13s  if  NT  is  even, 
NT/2+1  card  no.  13s  if  NT  is  odd. 


NOTE:  The  following  set  of  cards  specifies  the  catalytic  species  (Ml,  M2,  M3,  .  .  .  ) 
and  their  respective  composition  in  terms  of  the  species  participating  in 
the  reactions.  Weighting  factors  must  be  read  in  the  same  order  in 
which  the  thermodynamic  data  sets  are  read. 


14 

A6 

1-6 

AID(NS+1) 
(e.g..  Ml) 

Name  of  first  catalytic  species 

14a 

F5.2 

1-5 

WF(1,  1) 

Weighting  factor  of  first  species 
(for  first  catalytic  species).  Set 
weighting  factor  to  zero  for  any 
reactant  which  does  not  contribute 
to  the  respective  catalytic  species. 

F5.2 

6-10 

WF(1,  2) 

Weighting  factor  of  second  species 
contributing  to  first  catalytic  species. 

F5.2 

75-80 

WF(  1,  16) 

Weighting  factor  of  sixteenth  species 
contributing  to  first  catalytic  species. 

14a 

F5.2 

1-5 

WF(1,  17) 

Weighting  factor  of  seventeenth  species 
contributing  to  first  catalytic  species, 
etc. 

14 

A6 

1-6 

AID(NS+2) 

Name  of  second  catalytic  species. 

14a 

F5.2 

1-5 

WF(2,  1) 

Weighting  factor  of  first  species 
contributing  to  second  catalytic 
species,  etc. 

14 

A6 

1-6 

AID(NS  +NM) 

Name  of  last  catalytic  species,  etc. 

The  same  temperature  points  must  be  used  for  each  species. 


B-3 


NOTES:  1.  There  are  NM  card  no.  14s 

2.  There  are  as  many  card  no.  Has  as  needed  to  account  for  NS 
species  at  16  species/card 


Card  Format  Column  Variable 


Description 


A6 

1-6 

7 

A6 

8-13 

A6 

14 

6x 

15-20 

21 

A6 

22-27 

28 

A6 

29-34 

35 

A6 

36-41 

42-48 

12 

49-50 

11 

51 

E8.2 

52-59 

F5.2 

60-64 

F10.1 

54-74 

F6.2 

75-80 

8E10.6 

1-10 

PRINTX 

11-20 

X 

21-30 

DX 

31-40 

XMAX 

41-50 

XMIX 

Species  A 
+  sign 

Species  B  (or  M) 

+  sign 
Blank  or  M 
=  sign 
Species  C 
+  sign  (if  needed) 

Species  D  (or  M) 

+  sign  (if  needed) 

Species  E  (or  M) 

Blank 

Reaction  type,  1  to  12 
Rate  constant  type,  1  to  5 

A,  pre-exponential  factor  (cm- 
particle-sec)  units 

N,  temperature  exponent 

B,  activation  energy  (cal/mole) 
M,  temperature  exponent 

Output  print  increment  (cm) 
Initial  X  (usually  0.0  cm) 

Step  size  (cm) 

Terminal  station  (cm) 

Mixing  length  (cm) 


- Input  Cards  17,  18,  19  and  20  only  if  IGAS  =1 - 

NOTE:  Input  1  set  of  cards  17,  18  and  19  for  each  lasing  specie  (3  max) 


17 


8E10.6 


1-10 

WLM 

11-20 

WE 

21-30 

WEXE 

31-40 

BE 

41-50 

AE 

51-60 

RAS 

61-70 

RBS 

71-80 

SYMN 

—  Molecular  weight  of  the  lasing 
species 

System  constants  w  ,  u  x  ,  B  ,  and 
cte  (1/ cm)  of  the  lasing6 rnolecule 

Resonance  broadening  constant  a* 

(fC  2/ cm -atm) 

Resonance  broadening  constant  b* 
(K/ cm-atm) 

Symmetry  number 
=  1.0  for  diatomic  molecules 
=  2.0  for  CO 


^Listed  in  same  sequence  as  thermodynamic  species. 

Data  for  typical  lasing  molecules  may  be  found  in  Refs.  5  through  7. 
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Card 

F  ormat 

Column 

Variable 

18 

8E10.6 

1-10 

AB 

11-20 

BB 

21-30 

CB 

31-40 

AV 

41-50 

BV 

19 

2014 

1-4 

5-8 

JFIX(NV) 

etc 

Description 


Curve -fit  coefficients  for  the 
matrix  elements  of  the  dipole 
moment  for  the  v=l  —  v=0 
transition  of  the  lasing  molecule. 


JFIX(l)j  Lower  level  rotational 
JFIX(2)  quantum  numbers  for  all 
etc.  transitions  (1,NV)* 

If  input  as  zero,  the  program  will 
locally  select  all  J-values  based  on 
the  highest  gain  (J-shifting). 


NOTE:  The  following  card  group  specifies  the  collision  broadening  constants, 

aj(j  =  l»  NS),  one  for  each  species,  in  the  same  order  in  which  the  thermo¬ 
dynamic  data  sets  are  read  (input  only  if  LFLAG  =  1). 


20  8E10.6  1-10  RA(NS)  Collision  broadening  constant  for  each 

11-20  chemical  species, 

etc. 


NOTE:  Cards  21  and  21a  contain  the  initial  conditions  for  each  stream*  There 
is  one  set  of  cards  21  and  21a  for  each  stream 


21 


21a 


8E10.6 

1-10 

P 

Known  pressure  (torr) 

11-20 

U 

Known  velocity  (cm/sec) 

21-30 

T 

Known  temperature  (  K) 

31-40 

SIZ 

Known  stream  size  (cm) 

41-50 

51-60 

61-70 

71-80 

ALPHA(NS) 

Species  mole  fractions**  (species 

1  to  4) 

8E10.6 

1-10 

etc. 

ALPHA(NS) 

Seventh  species  mole  fraction,  etc. 

For  operation  at  fixed  J-values. 

**  • 

Listed  in  same  sequence  as  thermodynamic  tables. 
Stream  1  is  the  reacting/lasing  stream. 
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Card 

Format  Column 

Variable 

Description 

22 

15 

PPLOT 

Number  of  parameters  to  be 

plotted. 

23 

1615 

PLOT 

Input  control  flags  for  parameters 

to  be  plotted 
=  1,  power 
=  2,  pressure 
=3,  velocity 
=4,  temperature 
=  5,  density 

=xx6,  species  (where  xx  is  the 
particular  species) 

NOTE:  The  order  in  which  PLOT  is  input  is  arbitrary. 

NOTE:  Card  24  contains  the  coefficients  defining  the  width  of  the  lasing 

zone  and  is  read  only  if  ILASL  =  1.  The  width  is  given  by  a  second 
order  equation  of  the  form; 


ILASL  =  Qq  +  Q^x  +  Q^x^ 

24  3E10.6  1-10 

11-20 
21-30 


COEF(IO)  Q 
Q. 


0 
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Appendix  C 

ONE -DIMENSIONAL  LASER  AND  MIXING  PROGRAM 

(ODLAMP) 

(Sample  Cases) 

C-l  Single  Species  Lasing  (HF) 

C-2  Multi-Species  Lasing  (DF 

and  CO^) 


Appendix  C 

C-l  SAMPLE  INPUT 


UOLAMP  TEST  CASE  -  POwE*  ON  -  COMPARISON 
2  24  12.  SO  l  5  6  0 

20*0 
0.V5 
0.*b 
0.0 


0*0 


HF<0  ) 

20.008 

-65. 

,14 

o. 

6.967 

29.084 

-2.075 

100. 

6.950 

33.907 

-1 .379 

*00. 

6.961 

38.731 

-0*683 

300. 

6.964 

41 .554 

0.0  13 

bOO. 

6.970 

4b. 1 12 

1 .406 

700. 

6.97b 

47. 4b6 

2.611 

voo. 

6.986 

49.213 

4.197 

1200. 

7.007 

bl .225 

6.296 

1600. 

7.041 

b3.24b 

9.  10b 

2000. 

7.C80 

b4 . 82 1 

1 1.930 

2**  00  • 

7.  120 

b6. 1  lb 

14.770 

2800. 

7.160 

57.216 

17.626 

HF<  1) 

20.006 

-5  3. 

b  13 

0. 

6.05,7 

29.084 

-2.075 

1  00. 

6.958 

33. *07 

-1.379 

200. 

6.Q61 

38.731 

-0*663 

->0G. 

o.  9t>H 

41.054 

u*01j 

bOO. 

6.970 

4b. 1 12 

1  •  4Qo 

700  • 

6  •  9  7  O 

47 • 458 

2.61  1 

*00. 

6 . 9ttb 

49.213 

4.  197 

1200. 

7. 007 

bl  .22b 

6.296 

1600. 

7.041 

b3.24b 

9.  105 

2000. 

7 . 0  6  J 

54.621 

1 1 .930 

2400. 

7.120 

b6. 1  lb 

14.770 

2800. 

7.  160 

b7 . 2 1 6 

1 7.626 

HF  (  2  ) 

20.008 

-42. 

978 

0. 

6.957 

29.084 

-2.075 

100. 

6.95b 

90  7 

-1.379 

200. 

6.Q61 

38.731 

-0.666 

300. 

6 . 96  4 

4 1.584 

0.0  13 

500. 

6.970 

45. 1 12 

1 .406 

700. 

6.97b 

47.458 

2.611 

*00. 

6.980 

49.213 

4.  197 

1200. 

7. 0 07 

bl  .225 

6.296 

1600. 

7.041 

53.24b 

9.  105 

2000. 

7.080 

54.621 

1 1 .930 

2400. 

7.  120 

b6 .  1  lb 

14.770 

2800. 

7.160 

b7 .216 

17.626 

HF  (  J  ) 

20.008 

-32. 

622 

0. 

6.Q57 

29.084 

-2.075 

100. 

6. 95  o 

33.907 

-1.379 

200. 

6 . 96  i 

38.731 

-0 • 663 

300. 

6. 964 

4 1 . bb4 

0.0  13 

bOO. 

6.970 

4b. 1 12 

1  .406 

700. 

6.979 

47.458 

2.611 

voo. 

6.0  8  0 

49.213 

4.197 

1200. 

7.007 

bl .225 

6.296 

1600. 

7.041 

53.245 

9.  105 

2uG0  . 

7  •  o  b  J 

54.621 

1 1 . 9jO 

2400. 

7.12J 

56. 1 15 

14.770 

WITH  l-o  LAMP  AND  RESALE 

0  0  0  1  1  1  O  1 


50* 

6.957 

29*084 

-1 .727 

150* 

6.960 

36*726 

-1*031 

250* 

6.962 

40*284 

-0*335 

400. 

6.967 

43.557 

0.710 

600. 

6.973 

46 • 363 

2.  104 

800. 

6.961 

46 . 390 

3.499 

1000* 

6.992 

49.949 

4.696 

1400* 

7.023 

52.305 

7.699 

1800. 

7.061 

54.076 

10.516 

2200. 

7.  100 

55*497 

13.486 

2600. 

7.  140 

56.685 

16.196 

3000. 

7.  180 

57.71 1 

19.060 

50. 

6.957 

29 • 084 

-1 .727 

150. 

6.960 

36.726 

-1 .031 

250. 

6.962 

40*284 

-0.335 

400. 

6.98  7 

4j • 657 

0.710 

600. 

6.975 

46.363 

2.104 

600. 

6.961 

46.390 

5.499 

1000. 

6. 99£ 

49.949 

4.696 

1400. 

7.023 

52.300 

7.699 

1800. 

7.061 

54.076 

10.516 

*200. 

7.  100 

55.497 

13.486 

2600. 

7.  140 

56*685 

16.196 

3000. 

7.  180 

57.711 

19.060 

50. 

6.957 

29.084 

-1 .727 

150. 

6.960 

56 .726 

-1 *031 

250. 

6.962 

40*264 

-0*335 

400. 

6.967 

43*557 

0.710 

600. 

6. 973 

46 • 383 

2.104 

600. 

6.961 

46*390 

5.499 

1000. 

6.992 

49*949 

4.696 

1400. 

7.023 

52*305 

7  •  699 

1800. 

7.061 

54.076 

10.516 

2200. 

7.  100 

55.497 

1 3.466 

2600. 

7.  140 

56*685 

16.196 

3000. 

7.  180 

57.71  1 

19.060 

50. 

6.957 

29*064 

-1 .727 

1  50. 

6.960 

36.728 

-1.031 

250  . 

6.962 

40*284 

-0.335 

400* 

6.967 

43*557 

0*710 

600. 

6.973 

46 • 5bo 

2.104 

HOO. 

6.961 

48*390 

3.499 

1000. 

6.992 

49.949 

4.896 

1400* 

7.023 

52*305 

7.699 

1600. 

7.061 

54.076 

10.516 

c200. 

7.  100 

55.497 

15.468 

2600. 

7.  140 

56.665 

10.196 

C-l 


2600* 

7.160 

57.216 

17.62b 

3000* 

7.  160 

57*711 

19*060 

HF(-4) 

20 .006 

-22. 

7  32 

0. 

6.95  / 

29.064 

-2.075 

50* 

6.957 

29*084 

-1 *727 

1  00. 

6.958 

33*907 

-1 .379 

150* 

6.960 

36*728 

-1*031 

^00. 

6.961 

36.731 

-0*663 

250. 

6  •  962 

40*204 

-0 • o3b 

300. 

6.964 

4 1 .554 

0.0  13 

400* 

6.967 

43*557 

0*710 

500  • 

b.  970 

4b. 1 12 

1  *406 

600  • 

6.973 

4o  *383 

£•104 

700. 

6.976 

47.456 

2.611 

600* 

6.981 

46*390 

3.499 

voo. 

6.986 

49.213 

4.197 

1000* 

6.992 

49*949 

4  .  o90 

1200. 

7. 007 

51 .225 

6.296 

1400* 

7.023 

52*305 

7.699 

1600* 

7.041 

53.245 

9.  105 

1800* 

7.061 

54*076 

10*316 

2000. 

7.080 

54 . b2 1 

1 1.930 

2200. 

7.  100 

55*497 

10.406 

2400. 

7.  120 

56. 1 15 

14.770 

2600* 

7.  140 

56*663 

1 6 • 1 9o 

2b00  • 

7.  160 

57.216 

17.626 

3000* 

7.  180 

57*71 1 

19.060 

HF(b) 

20.00b 

-13. 

299 

o. 

6 •  Qp  7 

29.064 

-2.075 

50* 

6.957 

£9*064 

-1 .727 

1  00. 

6  •  9b  o 

-33.907 

-1.379 

150* 

6.960 

o6*72b 

-1.031 

200. 

6.961 

36.731 

-0*663 

250* 

6.96£ 

40*284 

-0.335 

3C0* 

6.964 

4 1 .554 

0*0  13 

400* 

6.967 

43*557 

0.710 

oOO. 

6.970 

45. 1 12 

1  .406 

600* 

6.973 

46 • 36o 

2.104 

700. 

6.97  b 

47.456 

2*81  1 

000  • 

6.901 

46*390 

o.499 

900* 

6.98b 

49.213 

4.197 

1000* 

6.992 

49*949 

4.096 

1200. 

7.007 

51 .225 

6.296 

1  400* 

7.023 

52*305 

7.699 

1600. 

7.041 

53.245 

9.  105 

1600* 

7.061 

54*076 

10*516 

£O00. 

7.080 

54 • o2 1 

1 1 . 9oO 

£200. 

7.  100 

33 .497 

1 O . 460 

2400* 

7.  120 

56. 1 15 

14.770 

2600* 

7.  140 

36 • 663 

10*196 

c'oOO  • 

7.  1  6  J 

57.216 

17.626 

3000* 

7.  1  BO 

57.71 1 

1 9 . G60 

HF(O) 

20.006 

-4  • 

313 

0. 

6.937 

29.064 

-2.075 

50. 

6.957 

29 • 0 6a 

-1 .727 

1  00. 

6.95b 

33.907 

-1 .379 

150* 

6.960 

36 • 72o 

-1*031 

200. 

6.^6  1 

36.731 

-0*663 

250* 

6.962 

40  •  £6*4 

— 0 • o3b 

JOC. 

6.964 

4  1  .  354 

0*0  13 

400* 

6.967 

4o . 557 

0*7  10 

500. 

6.970 

45. 1 12 

1  •  406 

600. 

6.973 

46 • 300 

£.104 

700. 

6  •  o76 

47.456 

2*611 

800* 

6.981 

46*390 

3.499 

-'GO  . 

b.Qdb 

49.213 

4.197 

1000* 

6*992 

49*949 

4.696 

1^00* 

7  .  CO  / 

b 1 . ^25 

6.296 

1400* 

7*023 

52  •  oO  -> 

7.699 

lbOO. 

7 . 04  1 

53 .24b 

9.  105 

1800* 

7*061 

54 • 07o 

16.516 

20  00  • 

7  .  o  8  0 

54.621 

1 1 .930 

'2200* 

7*  100 

33*497 

1 o . 466 

£^00* 

7. 120 

j6 .lib 

14.770 

2600* 

7*  140 

50*603 

16.1 96 

2  o  0  u  • 

7  .  l  6  J 

b7 . £  1 o 

17.626 

3000* 

7.  100 

57*711 

19.060 

f- 

IV. 

lo. 

34 

0* 

5  .  060 

29. 46  1 

-1.656 

50* 

5.066 

29*46 1 

-1 .495 

1  oo. 

b  .  C6  8 

32.  1  lb 

-1 .059 

150* 

5*325 

34  *20 1 

-0*003 

^oo. 

b  .  40  3 

35.746 

-0.534 

250* 

5*430 

36*960 

-0*262 

ooc . 

5.436 

37.95 1 

0*010 

400* 

5*361 

39*505 

0*550 

oOO. 

5.282 

40.693 

1  *082 

600* 

5*216 

4  1  .650 

1  .607 

/oo. 

=3.  1  O  V 

42.4b 

£  •  126 

600* 

5*  133 

43* 130 

£  *64  1 

J  0  0  • 

b.  103 

43.74 1 

3.  153 

1000* 

5*083 

44*277 

3*663 

1200  • 

5.032 

45.20 1 

4.676 

1400* 

5*032 

45*976 

5 . 684 

1  bOO  • 

5.018 

46.049 

6.669 

1600* 

5*009 

47*24 

7*692 

£000. 

b.  001 

47.767 

0.693 

£200. 

4.996 

46*244 

9.692 

2  <4  G  0  • 

4.992 

4©.o7o 

10*691 

2600. 

4.966 

49 • Q7o 

1 1 .009 

2oo0  • 

4 . 98b 

49.447 

12*687 

3000* 

4.964 

49*791 

13.683 

F2 

;vj. 

0. 

o. 

6.938 

35.67 1 

-2* 1 10 

50* 

6*958 

35*871 

-1 .762 

1  CO. 
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